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Abstract 

The equivalence of two optimality principles leading to Murray's law 
has been discussed. The first approach is based on minimization of bio- 
logical work needed for maintaining the blood flow through the vessels at 
required level. The second one is the principle of minimal drag and lumen 
volume. Characteristic features of these principles are considered. 

An alternative approach leading to Murray's law has been proposed. 
For that we model the microcirculatory bed in terms of delivering vascular 
network with symmetrical bifurcation nodes, embedded uniformly into 
the cellular tissue. It was shown that Murray's law can be regarded as a 
direct consequence of the organism capacity for controlling the blood flow 
redistribution over the microcirculatory beds. 

1 Introduction 

The great amount of natural systems have highly branching networks. As a 
evident example of such systems we may regard living tissue where blood sup- 
plies the cellular tissue with oxygen, nutritious products, etc. through branch- 
ing vascular network and at the same time withdraws products resulting from 
living activities of the cellular tissue. A similar situation takes place in respi- 
ratory systems where oxygen reaches small vessels (capillaries) going through 
the hierarchical system of bronchial tubes. There is a question of what physical 
principles govern network organization living systems under consideration. In 
this paper we focus our attention on the analyzes of different known optimal 
principles of network formation for microcirculatory bed and developing new 
approach to this problem. 
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Table 1: Typical parameters of the vessel arrangement of a 13-kg dog, after f^Of . 



Vessel type 


Diameter Length Tissue/ Vessel Arrangement 
2a, fim I, cm radii, d/a anisotropy* 


Primary arteries 
Small arteries 
Terminal vessels 
Arterioles 


300 1.0 30 0.86 
100 0.5 20 0.50 
50 0.2 10 0.37 
20 0.1 7 0.25 



■"measured as (nd 2 !) 1 ^ 3 jl 



2 Analysis of the optimality principles 

A microcirculatory bed can be reasonably regarded as a space-filling fractal 
being a natural structure for ensuring that all cells are serviced by capillaries 
[ pl| . The vessel network must branch so that every small group of cells, referred 
below to as "elementary tissue domain" , is supplied by at least one capillary. 
Since a typical length of capillaries is about 0.3 to 0.7 mm a vessel network 
generated by an artery of length of order of 1 to 5 cm should contain a sufficiently 
large number of hierarchy levels. At zeroth level we meet the host artery and the 
host vein, the mother and daughter vessels belong to n-th and (n + l)-th levels, 
respectively, and the last level N comprises capillaries. So at each level n of the 
vascular network the tissue domain supplied by a given microcirculatory bed as 
a whole can be approximated by the union of the tissue subdomains whose mean 
size is about the typical length l n of the n-th level vessels. Thus, the individual 
volume of these subdomains is estimated as V n ~ and their total number (as 
well as the total number of n-th level vessels) is about M n ~ Vo/l^, where Vb is 
the total volume of the microcirculatory bed. The higher is the level, the more 
accurate become the independence of such estimates from the particular details 
of vessel arrangements. For internal organs they approximately hold also for 
large vessels of regional circulation. To justify the latter statement we present 
Table |l| relating the vessel lengths and radii to the radii of the corresponding 
tissue cylinders, i.e. the cylindrical neighborhood falling on one vessel of a fixed 
level. 

This condition that the vascular network be volume-preserving from one 
generation to the next gives us immediately the local relation between the char- 
acteristic lengths of the vessels: I 3 w gl^+i (here g = 2 is the order of the 
vessel branching node). Whence it follows that l n ~ Z g - "/ 3 , where Iq is the 
characteristic size of the microcirculatory bed region or, what is practically the 
same, the length of the host artery. 

The following analysis, however, will require a more detailed information 
about the vascular network architectonics. Namely, we need to know how the 
vessel radii change at the nodes and the relative arrangement of mother and 
daughter vessels. Actually here we meet the problem as to what fundamental 
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Figure 1: Characteristics of the vessel branching 



regularities govern the vessel branching. These regularities manifest themselves 
in the relation between the radii a , a\, a 2 of mother and daughter arteries, 
respectively, and the angles 9\, #2, #12 = &i + #2 which the daughter branches 
make with the direction of the mother artery and with each other (Fig. [l]). 

A detailed theoretical attempt of understanding this regularity was made 
first by Cecil D. Murray in 1926 |l9, p0| . He proposed a model relating the 
artery radii at branching nodes (Fig. |) by the expression 
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with x = 3 (1) 



thereafter referred to as Murray's law (x is also called the bifurcation exponent). 
Then Murray's approach was under development in a large number of works, see, 
for example, @, |[ |§ |§ anda series of works by Zamir et all. M ^jj (40j [flf 



and by Woldenberg et all. [|33| (for a historical review see also 
The idea of Murray's model is reduced to the assumption that physiological 
vascular network, subject through evolution to natural selection, must have 
achieved an optimum arrangement corresponding to the least possible biological 
work needed for maintaining the blood flow through it at required level. This 
biological work V involves two terms: (i) the cost of overcoming viscous drag 
during blood motion through the vessels obeying Poiseuille's law, and (ii) the 
energy mctabolically required to maintain the volume of blood and the vessel 
tissue. Dealing with an individual artery of length I and radius a with a blood 
flow rate J in it we get: 

v= 8^! +m7ra 2 Z) (2 ) 
7ra 4 

where 77 is the blood viscosity and m is a metabolic coefficient. Minimizing 
function (|l|) with respect to a we find the relation between the blood flow rate J 
and the artery lumen radius a corresponding to the given optimality principles: 

J = ka 3 , (3) 
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Table 2: Integral characteristics of the vessel cross-section at different branching 



levels, after p). 



Vessel type 


Mean radius 
a, /im 


E« 2 Ea 3 -io 

cm 2 cm 3 


E« 4 

cm 4 








Homo sapiens 






Aorta 


12500 


1.56 1.95 


2.44 




Arteries 


2000 


6.36 1.27 


0.25 




Arterioles 


30 


127.4 0.382 


1.15 x 10 


-3 


Capillaries 


6 


1432 0.860 


5.16 x 10 


-4 


Venules 


20 


1273 2.55 


5.09 x 10 


-3 


Veins 


2500 


12.9 3.18 


0.80 




Vena cava 


15000 


2.25 3.38 

Canis familiaris 


5.06 




Aorta 


5000 


1.25 






Large arteries 


1500 


1.35 






Main arterial 
branches 


500 


0.75 






Terminal branches 


300 


0.486 






Arteriolas 


10 


0.400 






Capillaries 


4 


0.768 






Venules 


15 


2.70 






Terminal veins 


750 


7.59 






Main venous 
branches 


1200 


10.37 






Large veins 


3000 


10.80 






Vena cava 


6250 


1.53 
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where the coefficient k = yjrrm 2 / (I677) is a constant for the tissue under con- 
sideration. Due to the blood conservation at branching nodes we can write 
Jo = J\ + J2 (Fig. |l]) whence Murray's law (|l|) immediately follows. 

However, care must be taken in comparing measurements with prediction, 
particularly if averages over many successive levels are used. Already Mall 
himself noted that his data were approximate p7| . In particular, for large 
arteries of systemic circulation where blood flow can be turbulent the bifurcation 
exponent x should be equal to 7/3 « 2.33 as it follows from this optimality 
principle of minimum pumping power and lumen volume. 

There is also another optimality principle leading to Murray's law, the prin- 
ciple of minimum drag and lumen surface H?J . The drag against the blood 
motion through vessels is caused by the blood viscosity and can be described in 
terms of the shear stress on the walls of vessels, 2%alTt, where 

(4) 

for the laminar flow and 5P is the pressure drop along a vessel of length I 
and radius a. Then the given optimality principle is reduced to the minimum 
condition for the function 

where ml is a certain weighting coefficient. Minimizing (5) with respect to a we 



, leading to Murray's 
||, 0]) aimed at find- 



get a relationship between J and a of the same form as_(3) 
law again. There were a number works (see, e.g., E3, 
ing out what the specific optimality principle governs the artery branching by 
studying the angles of daughter vessels, 8\, O2, #12, in relation to the asymmetry 
of the branching node, 02/01 (Fig. |l|). However, on one hand, all the optimal- 
ity principles give numerically close relationships between the vessel angles and 
radii for the bifurcation exponent x ~ 3 |33fl . On the other hand, it turned out 
that experimentally determined branching angles generally exhibit considerable 
scatter around the theoretical optimum. The matter is that small variations of 
the total "cost" of artery bifurcation about several percents causes the actual 
vessel angles to deviate significantly from the predicted optimum. This feature 
is illustrated in Fig. ^ showing the variations in the vessel angles governed by the 
minimality of functional ([5]) with imposed 10% perturbations. Namely, varying 
the coordinates of the branching node (Fig. [I]) we get that the minimum of the 
function V'q -I- V'i + V' 2 is attained when 

a 1 cos 6*i + a 2 cos 9 2 = (1 + e)a , (6) 
01 sin 61 — a-i sin 62 = e'ao, (7) 

where the additional terms eao and e'ao with |e| , |e'| < 0.1 describes possible 
deviations from the optimality condition. The resulting values of 6\ and O2 
are depicted in Fig. ||. It should be noted that expressions (p|), (Q) correspond 
actually to the mechanical equilibrium of the node under the action of vessel 
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Figure 2: The angles that the daughter arteries make with the mother artery 
when the vessel branching is governed by the minimum drag & suface principle 
(solid lines) and under 10% perturbations (darkened region). 



walls strained by the blood motion and the additional terms describe a possible 
effect of the cellular tissue. 

Nevertheless, the optimality principles based on functional (Q) seems to gov- 
ern the artery bifurcations 0. Besides, this principles gives also adequate esti- 
mates of the integral characteristics of microcirculatory beds |L0| [ll| . 

The bifurcation exponent x, on the contrary, is well approximated by the 
Murray value, x « 3, at least starting from arteries of intermediate size J39| . 
This value meets also the space-filling requirement for the vascular network 
fractal in geometry to fill precisely the space of a fixed relative volume at each 
hierarchy level Q . Indeed, assuming the volume of the tissue cylinders match- 
ing an artery of length I and lumen radius a to be about Z 3 we get that the 
corresponding relative volume of blood is (a/l) 2 . So it is fixed if a = constant • I 
and, thus, <Zq = a 3 + a 3 , provided the tissue cylinder matching the mother artery 
is composed of the tissue cylinders of the daughter arteries. 

In order to specify the microcirculatory bed structure we need also to clas- 
sify vessels according to the symmetry of their branching (Fig. ||). The matter 
is that [|d| arteries with predominantly asymmetric bifurcations give off com- 
paratively little flow into its side branches along its way and, therefore, able 
to carry the mainstream flow across larger distances. Conversely, a more sym- 
metric bifurcation pattern splits flow into numerous small branches, thereby 
delivering blood to its surrounding tissue. Such arteries have been attributed a 
"conveying" and "delivering" types of function, respectively. Since blood must 
be conveyed towards the sites at which to be delivered, both types of vessels 
occur in real arterial trees. Moreover, a larger conveying vessel may switch into 
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a b levels 



Figure 3: Schematic illustration of the "conveying" (a) and "deliverying" (b) 
types of artery pattern. 

a bunch of small delivering branches. 

Obviously, real arterial trees should contain a great variety of intermediate 
stages in between these extremes and as our field of view moves from the large 
systemic arteries to small arteries of regional circulation the vessel bifurcation 
should become more and more symmetrical. This has been also justified by nu- 
merically modelling the structure of arterial trees governed by the minimality 
condition of blood volume . According to the experimental data 

(see, e.g., the work || and Fig. |4| based on it) even sufficiently large regional 
arteries of diameter and length about 300 /^m and 1 cm, respectively, (Table 
branch symmetrically, at least at first approximation.. Therefore microcircula- 
tory beds as they have been specified above can be regarded as a vessel network 
with approximately symmetrical bifurcations. 

In other words, we may think of the systemic arteries as vessels of the con- 
veying type where the mean blood pressure is practically constant. Conversely, 
the arteries of microcirculatory beds should belong to the delivering type and 
mainly determine the total resistance of the vascular network to blood flow, 
with the blood pressure drop being uniformly distributed over many arteries of 
different length. 



3 Physiological mechanisms governing the ves- 
sel arrangements 

The universality of Murray's law for distributed transport systems in many dif- 
ferent live organisms raises questions as to: What cues are available to organisms 
to use in generating such systems? What physiological mechanisms enable them 
to adapt to altering conditions? Do in fact live organisms follow certain global 
optimality principles? 

For Murray's law (^) the shear stress r t is constant (see formula (||)) through- 
out a given artery system. Rodbard proposed that the shear stress detected 
by the vessel endothelium leads to the vessel growth or contraction, and Zamir 
[ |38| suggested that this leads to the development of Murray's system as vessels 



7 




0.00 1 1 1 1 

10 100 1000 10000 
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Figure 4: Bifurcation symmetry vs. vessel diameter for the porcine coronary 
arterial tree (based on the data of [pi.) The darkness intensity indicates the 
density of the experimental points. 

maintain a constant value of shear stress. Concerning the particular mecha- 
nism by which organisms can implement the shear stress sensitivity we can say 
the following. Now it is established that the adaptation of conduit arteries as 
well as resistance arteries to acute changes in flow is mediated by the potent 
endogenous nitrovasodilator endothelium-derived relaxing factor, whose release 
from endothelial cells is enhanced by flow through the physical stimulation of 
shear stress (see, e.g., [Q and references therein). The adaptation of arterial di- 
ameters to long-term changes in the flow rate also occurs through a mechanism 
which appears to involve the sensitivity to shear stress and the participation of 
endothelial cells, but remains not to be understood well Q. 

It should be noted that the shear stress equality through a vascular network 
does not lead directly to a certain optimality principle. Different principles, 
for instance, (^J) and (^|), can give the same condition imposed on the shear 
stress. Moreover, it is quite possible that the case of this equality is of another 
nature. In particular, for large conduit arteries in the human pulmonary tree 
the bifurcation exponent x is reported to be in the range 1-2, whereas Mur- 
ray's law holds well beginning from intermediate conveying arteries [|4[ 
The matter is that in large systemic arteries the blood pressure exhibits sub- 
stantial oscillations because of the heart beating, giving rise to damped waves 
travelling through the systemic arteries. The value of the bifurcation exponent 
x = 2 matching the area-preserving law at the branching nodes ensures that the 
energy-carrying waves are not reflected back up the vessels at the nodes. How- 
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ever, this requirement is also can be derived from a certain optimality principle 

Summarizing the aforesaid we will model the microcirculatory bed in terms 
of a delivering vascular network with symmetrical bifurcation nodes embedded 
uniformly into the cellular tissue. Besides, the Murray's law will be assumed to 
hold. The latter is also essential from the standpoint of the tissue self- regulation, 
which will be discussed in detail in the next section. Here, nevertheless, we make 
several remarks concerning the given aspect too, because it could be treated 
as an alternative origin of Murray's law (^). Let us consider a symmetrical 
dichotomous vessel tree shown, e.g., in Fig. ||b. In order to govern blood flow 
redistribution over the microcirculatory bed finely enough so to supply with, for 
example, increased amount of blood only those regions where it is necessary and 
not to disturb other regions the blood pressure should be uniformly distributed 
over the microcirculatory bed, at least, approximately. 

The blood pressure drop SP n along an artery of level n (n = 0, 1, 2, ... , 
Fig. H) for laminar blood flow is 

w 

For the space-filling vascular network this artery supplies with blood a tissue 
region of volume about 1^ and, so, under normal conditions the blood flow rate 
J n in it should be equal to J n w jl^, where j is the blood perfusion rate (the 
volume of blood flowing through a tissue domain of unit volume per unit time) 
assumed to be the same at all the points of the given microcirculatory bed. 
Then formula (||) gives us the estimate 

ir \a, 

whence it follows that SP n will be approximately the same for all the levels, i.e. 
the blood pressure will be uniformly distributed over the arterial bed if the ratio 
l n /a n takes a certain fixed value, l n » constant-a„ and, thus, J„ rs constant'-a^. 
Due to the blood conservation at branching nodes we can write 

Jo - Ji + J 2 (9) 

(see Fig. |). The later gives us immediately Murray's law (|l|). In other words, 
Murray's law can be also regarded as a direct consequence of the organism capac- 
ity for controlling finely the blood flow redistribution over the microcirculatory 
beds. 



It should be noted that in the previous papers |14|, 15, [L6| we considered in 
detail the mathematical model for the vascular network response to variations 
in the tissue temperature on the given network architectonics. We have found 
that the distribution of the blood temperature over the venous bed aggregating 
the information of the cellular tissue state allows the living tissue to function 
properly. We showed that this property is one of the general basic characteristics 
of various natural hierarchical systems. These systems differ from each other 
by the specific realization of such a synergetic mechanism only. 
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